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Abstract 

We initiate a study of anisotropic plasma at strong coupling using 
the AdS/CFT correspondence. We construct an exact dual geometry 
which represents a static uniform but anisotropic system and find, that 
although it is singular, it allows for a notion of 'incoming' boundary 
conditions. We study small fluctuations around this background and 
find that the dispersion relation depends crucially on the direction of 
the wave- vector relative to the shape of the anisotropy reminiscent of 
similar behaviour at weak coupling. We do not find explicit instabili- 
ties to the considered order but only a huge difference in the damping 
behaviour. 



1 Introduction 

One of the outstanding problems in our understanding of heavy-ion colli- 
sions is the question why quark-gluon plasma (QGP) can be described very 
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accurately by hydrodynamics so soon after the collision pQ . The difficulty in 
addressing this question lies partly in a mixed weak/strong coupling physics 
of the initial state and early dynamics. Although some notable work has 
been done (see e.g. [2]), we are still far from a definite understanding. 

Let us recall that hydrodynamics by definition involves the notions of 
(isotropic) pressure while just after the collision the energy-momentum tensor 
is definitely anisotropic. In a first approximation the question is thus to 
understand the process of isotropisation of an anisotropic plasma system. 

One of the mechanisms which was suggested to be responsible for this 
behavior is the appearance of instabilities in an anisotropic plasma system 
first discovered in [3J. Since the study of an expanding (anisotropic) plasma 
system is quite complicated if not impossible, a simpler system has been 
investigated namely an anisotropic plasma sytem which fills the whole space 
and evolves in Minkowski time (and not proper time). 

The instablities at weak coupling have been investigated in detail [U |5| . 
Subsequently the process of isotropisation has been studied in real time 
through numerical simulations [U El 0, E] • Initially, for weak fields the evo- 
lution is exponential, in accordance with the instabilities discovered in weak 
coupling computations around an anisotropic plasma background, then when 
nonlinear field effects become stronger, the evolution becomes linear in time. 

The motivation of this paper is to address similar issues at strong coupling 
using the AdS/CFT correspondence [H] as a calculational tool. Our ultimate 
goal is to study the temporal evolution of the anisotropic plasma system 
and its approach to isotropy. However in the present paper we want to 
first investigate the situation when the anisotropic plasma is assumed to be 
static and to look for possible instabilities of small fluctuations in direct 
analogy to the weak coupling considerations. We plan to study the real-time 
isotropisation process in future work [10J. 

The plan of this paper is as follows. In section 2, we will briefly review 
some features of plasma instabilites at weak coupling. In section 3 we review 
the AdS/CFT framework used and in section 4 we construct the geometry 
dual to a static anisotropic plasma system. In the following two sections 
we comment on some of its pathologies and discuss the issue of defining 
physically natural boundary conditions. In section 7 we study the disper- 
sion relation of R-charge fluctuation modes and we close the paper with a 
discussion and an appendix containing the relevant wave equations. 
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2 Plasma instabilities at weak coupling 



In this section we will briefly describe some qualitative features of plasma 
instabilities studied at weak coupling. The majority of numerical and analyt- 
ical work has been done for an infinite, spatially uniform system which does 
not expand. One basically starts from some initial conditions when momen- 
tum distributions are anisotropic. Typically one has a separation between 
hard and soft modes and studies how the hard modes lead to isotropisation 
of the soft modes but this is not strictly necessary [8]. Then one studies 
the time dependence of electric/magnetic fields and consequently the time 
evolution of the energy-momentum tensor: 



T 



( e \ 

p L (t) 

pr(t) 

\ pr(t) J 



(1) 



where e = pl + 2pr- Such studies have to involve numerical computations. 

Alternatively, one can try to compute the poles of the gluon propagator 
in the anisotropic system with a momentum distribution 

f(p) = yfiTt fi~(? + &i) (2) 

where fi SO is some isotropic distribution and £ is the anisotropy parameter 
related to the ratio of transverse to longitudinal pressure through 

t = - ~ 1 (3) 

PL 

The outcome is that some modes develop unstable behavior. The precise 
behavior crucially depends on the sign of the anisotropy. If £ > 0, then 
all modes with transverse wave vectors remain stable, while the longitudinal 
ones develop an instability in a finite range of kj,. If £ < the situation 
is reversed with the longitudinal modes remaining stable and the transverse 
modes developing an instability. 

The unstable modes can be roughly identified with the initial behavior 
of the plasma in the time-dependent simulations outlined above. Thus the 
simpler computation of the modes gives an indication of the direction of the 
evolution of the system. Later when nonlinear effects become important, the 
evolution ceases to be exponential. 
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In the strong coupling regime one may attempt to study these questions 
in both ways. In this paper we will analyze what happens when we try to 
mimick the simpler approach, namely to consider an anisotropic system and 
look at small fluctuations. We plan to analyze the more realistic case of time 
dependence in the future. 



3 The AdS/CFT framework 

Within the AdS/CFT correspondence, a system of plasma is described by 
a dual geometry which is constructed as follows. Let us assume that the 
only nonvanishing expectation value in the system is a certain profile of the 
energy-momentum tensor: 

OVM> (4) 

which has a specific dependence on the Minkowski coordinates x M . Then the 
dual geometry is written in the Fefferman-Graham coordinates as 

_ 2 Q fj J i/['^^ ' i z^jdoc^ doc ~\~ dz> fQ"^ 

z 2 

The metric is determined by solving Einstein's equations with negative cos- 
mological constant (A = —6) which can be written as 

R*p + ^g a p = o (6) 

with the boundary condition at z = given by 

9fiu (x^ z) = n,., + z 4 g$(x») + O (/) (7) 

where g$(x^) is related to the expectation value of the gauge theoretical 
energy-momentum tensor [11J: 

N 2 

(T^M) = 2^£tV) (8) 



The proposal put forward in p.2j is to determine the geometry for a given 
energy-momentum profile and choose the physical one by requiring the non- 
singularity of the resulting bulk geometry. 

This framework is geared to study the time-dependent processes and thus 
can be used to study isotropisation in real time starting from the time- 
dependent energy-momentum tensor ([T]) by constructing the dual geometry 
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and requiring nonsingularity. We plan to follow this route in future work 

ma. 

Because the above procedure, although conceptually simple, is technically 
quite involved, we decided to analyze what are the fluctuations around a 
system where we neglect the time dependence of the anisotropy of the plasma. 
This problem is a direct analog of the computation of the (unstable) poles 
of the gluon propagator in the unstable medium. Similarly as in the weak 
coupling case we do not expect such a static system to exist indefinitely. 
Certainly even if the system would be classically fine-tuned to stay in the 
unstable regime, quantum fluctuations would cause it to evolve. Therefore 
we expect the geometry to be somewhat pathological. The question that we 
wanted to ask is whether any kind of information may be extracted from it 
and see whether the pattern of fluctuation modes has some resemblance to 
the weak coupling situation. 



4 Geometry dual to a static anisotropic plasma 
system 

Let us now find the dual geometry corresponding to a uniform, static and 
anisotropic energy- momentum tensor: 



( e \ 

Pl 

p T 

\ p T ) 



(9) 



with e = pl + 2pt- The most general metric with these symmetries has the 
form 

(10) 



ds 2 = — (-a(z)dt 2 + b(z)dx 2 L + c{z)dx\ + dz 2 ) 



The scalar functions a(z), b(z) and c(z) have to vanish at z — 0, and their 
z 4 coefficients are related to the transverse and longitudinal pressure. The 
general solution of (jUJ) satisfying these constraints reads 



a{z) 
b(z) 
c(z) 



'I + A 2 z i )^~* VMZ2m (l - A 2 z i )^ v ^ l2W 



1 + A z 



1 - A Z 
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(11) 

(12) 
(13) 
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and the A and B parameters are related to the energy density and pressure 
through 

e = ^A 2 VS6 - B 2 (14) 



p L = -AV36 -B 2 - -A 2 B (15) 
o o 

p T = Ia 2 V36 - B 2 + l -A 2 B (16) 

6 3 

It is also convenient to link the B parameter with the £ anisotropy parameter 
defined through ([3]) to get 

B = ^ (17) 

v/18£ 2 + 48£ + 36 

When there is no anisotropy, B = and the above solution reduces to the 
standard static AdS black hole solution. 



5 The issue of the singularity 

An unavoidable property of the metric is, that once we have nonzero anisotropy, 
there is a singularity in the bulk, which is a very significant obstacle towards 
making a physical interpretation of our setup. This can be interpreted that 
a static anisotropic plasma system cannot exist at strong coupling. In fact 
such a statement is also true at weak coupling as discussed in section 2. 

We would like to try to interpret our setup as a snapshot of the dynami- 
cal evolving scenario (close to the initial condition) and to investigate small 
fluctuations around such a system in order, eventually, to compare with a 
similar setup at weak coupling where the pattern of stable/unstable modes 
has a specific dependence on the sign of anisotropy and the relative orienta- 
tion of the wave- vectors. 

This might give a hint towards the real time-dependent evolution of the 
anisotropic system and what would be the differences in behavior with what 
is known at weak coupling before attempting a real dynamical computation 
at strong coupling. 

Let us note, however, that because of the appearance of the singularity 
it is not clear what is the regime of validity of these results. The 'snapshot' 
interpretation could perhaps be validated by assuming very big 'classical' 
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occupation numbers but we do not know how to estimate possible timescales 
appearing at strong coupling. We decided to proceed nevertheless with the 
computation of small fluctuations and to see a-posteriori whether the ob- 
served behavior is physically viable, or whether it is obviously pathological. 

6 Boundary conditions 

One of the major conceptual issues that one has to face when considering 
spacetimes with naked singularities is the problem of what boundary con- 
ditions to impose at the singularity. This problem has been considered in 
the general relativity literature [T5| [H] , and not surprisingly the problem of 
stability/instability of such spacetimes like the negative mass Schwarzschild 
black holes depends on the chosen boundary conditions. 

In fact even for the ordinary Schwarzschild geometry if one would choose 
outgoing boundary conditions at the horizon we would get unstable modes. 
This issue of course never arises in the case of the black hole horizon, as there 
incoming boundary conditions are clearly singled out physically. 

We will try to find an analog of the incoming boundary conditions for 
our geometry and use it for our computations. In fact the geometry is not so 
pathological as e.g. a negative mass black hole where a similar construction 
would not be possible. This feature of ffTUj) is somewhat reassuring. 

Let us consider the wave equation for a massless scalar field in the ge- 
ometry (flQj) . For simplicity we will set A = 1. Performing the standard 
separation of variables 




.3 




and a subsequent change of variables 



x = 



- arctanh z 4 
4 



(19) 



we obtain the scalar equation in the form 




(20) 



7 



where the singularity is at x = oo. We see that close to the singularity the 
piece proportional to uo 2 dominates and the equation looks there as 



d 2 <f) 



+ 8u; V 2 ^-^"^ = 



(21) 



dx 2 



For zero anisotropy B = this has the familiar e i ^ ux ) e +iV8ux so i u tj ons 



which correspond to incoming and outgoing waves respectively. If we would 
consider the limit of very small B, we see that the solution behaves as for an 
ordinary horizon, with a clear separation of incoming and outgoing waves, 
which would only be modified very close to the singularity. In fact linearizing 
the equation in B will reduce the boundary conditions to the black hole ones. 
For general B, the solution of (12ip is a combination of Hankel functions 



with C = 6 — V36 — 2B 2 . As is well known Hankel functions form a conve- 
nient basis of incoming and outgoing wavefunctions in a cylindrical geometry. 
Hence in general we also have a clear definition of incoming and outgoing 
waves. 

Since we do not want any information flow from the singularity into the 
bulk, we will always pick these 'generalized' incoming waves as our boundary 
condition at the singularity. 

Let us note that the geometry ( jTUl) is rather special in that it allows for 
such a choice to exist at all. If we were to perform the same analysis, say for a 
negative mass (planar) AdS black hole we would have found that for nonzero 
wave-vectors these terms would dominate the u 2 term and give solutions for 
which there would be no notion of an incoming/outgoing wave interpretation. 

7 R-charge fluctuation modes 

In view of the appearance of plasma instabilities at weak coupling we are 
interested mainly in studying low lying fluctuation modes of the geometry 
(jTUl) . Let us therefore briefly review the situation for isotropic plasma [J5j 
[T6l [T7j . There all the modes exhibit damping (quasinormal modes). The 
smallest damping is associated with modes in the hydrodynamical regime. 




(22) 
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These involve modes of the graviton in the shear and sound channels which 
have the dispersion relations 



n ,, 1 , 2 k 2 
uj = -% — — k 2 + ... uu = —k-i-ri h... (23 

e+p ^3 3 e + p 

respectively, as well as modes of the bulk gauge field associated with corre- 
lation functions of R-charge currents 

u = -iD R k 2 + ... (24) 

where D R = 1/2ttT is the R-charge diffusion constant. All other modes have 
finite non-zero damping {Imu < const < 0) even for very long wavelength 

modes [mmE]. 

If we are to see an instability we have therefore to concentrate on the 
modes for which the (negative) imaginary part is smallest i.e. on the hydro- 
dynamic modes. In the present paper we will study the bulk gauge fluctua- 
tions as these are technically simplest. In the major part of the paper we will 
set A = 1 and reinstate it in the final results and formulas in the appendices. 
In these units the dispersion relation for R-charge diffusive modes (124"I) for 
isotropic plasma i.e. with B = takes the form 

" = - i 2^ + - (28) 

In this section we will study how this dispersion relation is modified by the 
presence of anisotropy. 

The equation of motion for the bulk gauge field takes the form 

8 a (V^F a ?) = (26) 

Due to the anisotropy of the plasma configuration and hence of the bulk 
geometry ( flOl) . one of the spatial coordinates is singled out (the longitudinal 
one), we will denote it by y, while the two transverse ones will be denoted by 
Xi t z- Similarly we will have a decomposition for wave vectors in the Fourier 
transform of the vector potential, 

A ^ X) = J ^ e ^ t+CkSA ^k) (27) 
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We will separately study the extreme situations when the wave vector is 
purely longitudinal or purely transverse: 



L : fc = (fc L ,0,0), T : q = (0, 0, k T ) (28) 

Let us recall that at weak coupling only one type of these modes develops an 
instability depending on the sign of the anisotropy parameter. 

Furthermore, since we are dealing with vector fields so there will be modes 
of the fields parallel and orthogonal to the wave vector. As a result we get 
two sets of coupled differential equations which can be further simplified by 
introducing gauge invariant variables, different in each set: 
Longitudinal modes: 

(L-L) E y (k L , z) = uA y (k L , z) + k L A t (k L , z), k L \\E y (29) 
(L-T) E 1 (k L , z) = uA^kL, z), E 2 (k L , z) = ujA 2 (k L , z) (30) 

Transverse modes: 



(T-T) E^kr, z) = uA 1 (k T , z) + k T A t (k, z), hr\\E u (31) 
(T-L) E y (k T , z) = uA y (k T , z), E 2 (k T , z) = ooA 2 (k T , z) (32) 

After that we are left with five quite lengthy equations. In the isotropic limit 
of the standard static black hole it is the electric modes which are parallel 
to the wave- vector which exhibit diffusive behavior (here (T-T) and (L-L)). 
The rest have O (1) damping even for small k. 

In the present paper we will be mostly interested in effects appearing for 
small anisotropy, hence we will linearize the resulting equations in B. We 
quote all resulting equations in this regime in Appendix A. 

Here we will analyze the equations for (T-T) and (L-L) setting A — 1, 
reinstating A dependence in the final answer. 

Close to the singularity z — 1, in the approximation linearized in B 
we find the following asymptotic behaviors which can be associated with 
incoming boundary conditions 

E y ~ (1 - zy' l ^ + ^ E x ~ (1 - z)'^'^ (33) 

Let us now concentrate on the equation for E y and make a decomposition 

Ey(z) = (1 - zr& + *g(u) (34) 
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where we used the variable u = z 2 . Furthermore, as we are interested in 
small frequencies and wave vectors we will rescale u — > euo and ki — *> ek^ 
and write 

g{u) = 1 + eg a (u) + e 2 g b (u) + B{£{u) + eg\{u) + ...) + ... (35) 

We then impose vanishing conditions at u — 1 for the gf(u)'s which deter- 
mines them completely. At this stage the mode has the correct incoming 
boundary condition at the horizon. Imposing Dirichlet conditions at the 
boundary gives the dispersion relation for the modes: 

g(0) = 1 + g a (0) + g b (0) + B(#(0) + <tf(0)) = (36) 

We find for the leading terms 



and for the subleading ones 



(rf(O) = ^(7r 2 -121og 2 2V (40) 



Let us analyze the condition (I3"6"j) . For zero 5, we recover immediately 
from (I37p . Once we turn on a small anisotropy, we find that for small k^ we 
cannot sustain the scaling u oc k 2 . Indeed then the first term in (I38I) will 
start to dominate and we get essentially 

1-^ = (41) 

Here we see that the outcome depends crucially on the sign of the anisotropy. 
If B > then, although we get no instability, the damping vanishes to this 

order and one has 

[~B 

uj = \ —k L + ... (42) 
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On the other hand if B < we get very strong damping^ 



U = ~ % \l~f kL + --- (43) 

When we increase ki the first term in ( l37j) -()40l) that will start to become 
important will be the first term in (1371) when k^ becomes of of order y>B|. 
At this stage we have to take into account both terms and we get 

h 2 ih 2 

1-^5 + ^ = (44) 

This is a quadratic equation which can be solved to give (after inserting 
dependence on A which was previously set to 1) 



«, = i V V ~ ( 45 ) 

We see that for B < 0, the sign of the square root has to be chosen as above 
in order to reproduce smoothly the B — > limit. This justifies the choice of 
sign in (|43|) . For positive B of course both signs are possible and there is a 
branch point on the real axis. 

Let us note that the change in the qualitative behavior depending on the 
sign of the anisotropy is quite similar to the one seen at weak coupling. For 
B > 0, at weak coupling modes with longitudinal wave vectors develop an 
instability while for B < they remain stable. Here the longitudinal ones 
develop a linear regime for B > while for B < they exhibit very strong 
damping. 

For modes with nonzero transverse modes the situation is reversed. The 
equation that gives the dispersion relation is: 

k 2 ik 2 

1 + 7^ + 4^ = 0, (46) 
and the solution reads (after reinstating the dependence on A): 



k 2 T + v k\ + \BA\k\ 
AV2A^ 



1 The sign is chosen as to reduce to the isotropic case when taking more terms into 
account and performing B — * 0. 
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In this case the term with B differs by a numeric constant but, what is more 
important, also has a different sign. This means that the behaviour of the 
modes with transverse wave vectors has an opposite pattern of strong/weak 
damping relative to the sign of the anisotropy. This is qualitatively similar 
to the weak coupling behaviour reviewed in section 2. 

One can also compute two-point functions using the standard procedure 
and find that the poles coincide with the dispersion relations found above. 

8 Discussion 

In the present paper we have studied the properties of a uniform infinite 
plasma system with anisotropic pressure. We have constructed the corre- 
sponding dual geometry and found that it is singular. This shows that such 
a system cannot be considered to exist indefinitely. The singularity is rela- 
tively mild in the sense that for a scalar wave equation in this background 
geometry there is a natural notion of ingoing boundary conditions in contrast 
to generic singular geometries. This feature allowed us to investigate small 
fluctuations and to look for possible instabilities and whether they would 
manifest themselves at the linear level in the supergravity description. 

Similarly as for weak coupling, we find that the qualitative behaviour of 
longitudinal and transverse modes depends crucially on the sign of anisotropy. 
In contrast to weak coupling however, to leading order in the anisotropy pa- 
rameter, we do not find instabilities but massless-like propagation (albeit 
with some quadratic damping) for modes with wave vectors for which one 
would expect instabilities at weak coupling, and very strong damping for 
modes which would remain stable at weak coupling. This qualitative cor- 
relation between the sign of the anisotropy and behaviour of the modes is 
reassuring in view of our doubts on using the singular geometry at all. Per- 
haps the lack of instabilities seen at the linearized level might be analogous to 
the lack of exponential growth for strong fields in the numerical simulations 
at weak coupling. In addition, Boltzmann-Vlasov simulations with large col- 
lision rate (corresponding to strong coupling) also showed the disappearance 
of instabilities |20j. 

One could extend the present investigation by going to higher orders in 
the anisotropy parameter but the most interesting question in our opinion 
is to determine the real-time evolution of the anisotropic system at strong 
coupling. This dynamical problem, although technically much harder, would 
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not be plagued by conceptual issues dealing with singularities. It would 
be very interesting to confront the outcome with the linearized results of 
the present paper. We plan to investigate the real-time evolution in future 
work [TO] . 
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A Linearized equations for R-charge modes 

In the following equations we neglected quadratic terms in B. 
Longitudinal wave-vectors 

3 + 4y/Z(3 + B)z 4 + 9Az 8 
El 3z(l-AzS) E * (48) 



+ - *2U ~ v^ 4 r% + VAzr^E, = 0, (49) 



3fc£(l - VAz*) 2 (l + VAz*)f(l - 3VAz\A - Viz*)) 
A, : — — 1= — : — : 1= — ~s : 1= — 1=. — : r 



3z(l - Az 8 )(k 2 L {l - VAz 4 ) 2 (l + VAz A p - u 2 (l - y/Az*)%(l + ^Az 4 ) 2 ) 

(50) 

u 2 {l - a/Iz 4 )! (1 + \/Az 4 ) 2 (3 + \[Az\Yl -8B + 9^/Az 4 )) 

32(1 - Az s )(k 2 (l - y/Az*) 2 (l + VJz 4 )f - u 2 (l - y/Az*)f(l + y/Jz 4 ) 2 ) 

(51) 

+ ( ^ (1+ ^ 4) - kl(l - VAzTHl + VAzr i+ ^)E y = (52) 

v (i - Vaz 4 ) 2 m ' y y ' 

(equation for E 2 is exactly the same as for the E 2 ). 
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Transverse wave-vectors 

3 + 4^(3 + ^ + 9^ 
2 3z(l-Az 8 ) 2 K ' 

+ - ^ - ^ 4 )-f (1 + VI*)- 1 -* )E 2 = 0, (54) 

(l — \j Az^y 



3 + v^ 4 (12-8S + 9v / As 4 ) F/ 
^ 3z(l - li^ ^ (55) 



+ ( _ JLJ - 4(1 - >/A8 4 )-t(1 + VAz*)- 1 --^ = 0, (56) 



E n , -3^(1 - v^) 2+ *(l ~ 3y^ 4 (4 - v^, 4 )) 
1 3^(1 -Az 8 )(k*(l - y/Iz*) 2+ T -uj 2 {1 - VAz*) 2+ ^ 



e 



u\l + v^ 4 ) 2+f (3 + 4^(3 + B)z* + 9Az s ) 



3*(1 - Az 8 ){k 2 T (1 - v 7 ^ 4 ) 2 ^ - cj 2 (1 - V^z 4 ) 2 "^)' 

^ 2 (l + v^ 4 ) _ ^ _ ^ 4) -f (1 + v^/)-i-f = o (59) 
V (1 - v 7 !^) 2 
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